Observation of quantum state collapse and revival due to the single-photon Kerr effect 
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Photons are ideal carriers for quantum information as they can have a long coherence time and can 
be transmitted over long distances. These properties are a consequence of their weak interactions 
within a nearly linear medium. To create and manipulate nonclassical states of light, however, 
one requires a strong, nonlinear interaction at the single photon level. One approach to generate 
suitable interactions is to couple photons to atoms, as in the strong coupling regime of cavity QED 
systems [1, 2]. In these systems, however, one only indirectly controls the quantum state of the light 
by manipulating the atoms [3]. A direct photon-photon interaction occurs in so-called Kerr media, 
which typically induce only weak nonlinearity at the cost of significant loss. So far, it has not been 
possible to reach the single-photon Kerr regime, where the interaction strength between individual 
photons exceeds the loss rate. Here, using a 3D circuit QED architecture [4], we engineer an 
artificial Kerr medium which enters this regime and allows the observation of new quantum effects. 
We realize a Gedankenexperiment proposed by Yurke and Stoler [5], in which the collapse and 
revival of a coherent state can be observed. This time evolution is a consequence of the quantization 
of the light field in the cavity and the nonlinear interaction between individual photons. During 
this evolution non-classical superpositions of coherent states, i.e. multi-component Schrodinger cat 
states, are formed. We visualize this evolution by measuring the Husimi Q-function and confirm 
the non-classical properties of these transient states by Wigner tomography. The ability to create 
and manipulate superpositions of coherent states in such a high quality factor photon mode opens 
perspectives for combining the physics of continuous variables [6] with superconducting circuits. 
The single-photon Kerr effect could be employed in QND measurement of photons [7] , single photon 
generation [8], autonomous quantum feedback schemes [9] and quantum logic operations [10]. 



A material whose refractive index depends on the inten- 
sity of the light field is called a Kerr medium. A light 
beam traveling through such a material acquires a phase 
shift ^Kerr — Kt I [11] whcrc / is the intensity of the beam, 
T is the interaction time of the light field with the ma- 
terial, and K is the Kerr constant. The Kerr effect is 
a widely used phenomenon in nonlinear quantum optics 
and has been successfully employed to generate quadrature 
and amplitude squeezed states [12], parametrically convert 
frequencies [13], and create ultra-fast pulses [14]. In the 
field of quantum optics with microwave circuits, the direct 
analog of the Kerr eflFect is naturally created by the non- 
linear inductance of a Josephson junction (specifically the 
0^ ^ (6 + b^)^ term in the Taylor expansion of the coscj) 
of the Josephson energy relation) [15, 16]. This eflFect has 
been used to create Josephson parametric amplifiers [17-19] 
and to generate squeezing of microwave fields [20]. How- 
ever, in both the microwave and optical domains, most ex- 
periments utilize the Kerr nonlinearity in a semi-classical 
regime, where the quantization of the light field does not 
play a crucial role. The Kerr eflFect for a quantized mode of 
light with frequency c^c, can be described by the normal or- 
dered Hamiltonian i^Kerr = hhJcCi^ a — h^a^ o) aa with K the 
Kerr shift per photon [2, 16]. The average phase shift per 
photon is again given by ^Kerr = K/ n with n the decay rate 
of the photon mode. Typical Kerr effects are so small that 
they are not visible on the single photon level as > 
Applications which require K much bigger than n include 
the realization of quantum logic operations [10], schemes 
for continuous variable quantum information protocols [6] 



and quantum non-demolition measurements of propagating 
photons [7]. 

The Kerr nonlinearity of a Josephson junction is also 
routinely employed to create superconducting qubits. In 
this case, one engineers a circuit with such a strong an- 
harmonicity that it can be considered as a two level sys- 
tem. By combining a qubit with linear resonators one re- 
alizes the analogue of strong coupling cavity QED, known 
as circuit QED. [1] This can be used to protect the qubit 
from spontaneous emission, manipulate and readout its 
quantum state and couple it to other qubits. One con- 
sequence of coupling any resonator to a qubit, is that the 
resonator always acquires a finite anharmonicity be- 
coming a Kerr medium itself [15]. In this paper, we have 
designed K large enough {K/ k > 30) to be well within the 
single-photon Kerr regime. At the same time, the non- 
linearity is small enough that short displacement pulses 
(^puise ~ 10 ns << 1/K) applied to the resonator create a 
coherent state, allowing us to conveniently access the large 
Hilbert space of the oscillator. Only a few experiments have 
previously come close to the limit where K/k ^ 1 while still 
maintaining the ability to create coherent states [21-23]. 

As a first demonstration within this single-photon Kerr 
regime, we realize a proposal [5] for creating photonic 
Schrodinger cat states. Specifically, we generate coherent 
states with a mean photon number of up to four photons 
and measure the Husimi Q-function of the resonator state 
using a new experimental measurement protocol. We then 
show the high quality of the Kerr resonator by measur- 
ing the time evolution of the collapse and revival of a co- 
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FIG. 1. Device layout and energy level diagram of the 
two cavity, one qubit device, a, Photograph of one half 
of two aluminum 6061 waveguide cavities coupled to a verti- 
cal transmon qubit. The right-hand side of the figure shows a 
detail of the qubit, fabricated on a c-plane sapphire substrate 
1.4 mm wide, 15 mm long, and 430 microns thick. The coupling 
strength of the qubit is determined by the length of the stripline 
coupling antenna which extends into each cavity. The upper 
cavity, with a resonant frequency of cJm/STT = 8.2564 GHz, is 
used for qubit readout and the lower cavity, with a frequency of 
uJc/2ti — 9.2747 GHz, is used to store and manipulate quantum 
states, b. Combined energy level diagram of the qubit coupled 
to the storage cavity. The qubit states are denoted as \g) and 
|e) respectively while the cavity states are labeled as |n) with n 
the number of photons in the cavity. Each photon in the cav- 
ity reduces the qubit transition frequency by x- Equivalent ly, 
exciting the qubit reduces the cavity transition frequency by x- 
The energy levels of the cavity are not evenly spaced due to the 
induced Kerr anharmonicity K/2n — 325 kHz. 

herent state. During the evolution, highly non-classical 
superpositions of coherent states, i.e. multi-component 
Schrodinger cat states, are formed, which show the coherent 
nature of the effect. This revival, in contrast to the Jaynes- 
Cummings revival of Rabi oscillations of a qubit induced by 
a coherent state [2, 24], is the revival of a coherent state in 
a resonator. An analogous effect was indirectly observed in 
early experiments with a condensate of bosonic atoms in an 
optical lattice [25]. Additionally, we confirm the nonclassi- 
cal properties of the transient states by performing Wigner 
tomography. 

We experimentally realize a highly coherent Kerr 
medium by coupling a superconducting "vertical" trans- 
mon qubit to two 3D waveguide cavities as shown in Fig. la. 
This design is based on a recently developed 3D circuit 
QED architecture [4]. The two halves of the cavities are 
machined out of a block of superconducting aluminium 
(alloy 6061 T6). Both cavities have a total quality fac- 
tor of about 1 million, corresponding to a single photon 
decay rate = 10 kHz. The vertical transmon con- 

sists of a single Josephson junction embedded in a trans- 
mission line structure which couples the junction to both 



cavities. The observed transition frequency of the qubit is 
ujq/27T = 7.8503 GHz and its anharmonicity is Kq/27r = 
(cjge — co'ef)/27r = 73.4 MHz using the standard convention 
for labeling from lowest to highest energy level in the qubit 
as (g,e,f,h,...). The energy relaxation time of the qubit is 
Ti = 10 /is with a Ramsey time T2* = 8 /is. The qubit 
is used to interrogate the state of the storage cavity which 
acts as a Kerr medium. The other cavity is used to read 
out the state of the qubit after the interrogation, similar to 
ref. [26]. 

The analysis of the distributed stripline elements and 
the cavity electrodynamics can be performed using finite- 
element calculations for the actual geometry. Combined 
with "Black-Box" circuit quantization, [15] one can derive 
dressed frequencies, couplings, and anharmonicities with 
good relative accuracy (see supplementary material). For 
the purposes of the experiments discussed here, the cou- 
pling of the qubit to the storage resonator, in the strong 
dispersive limit of circuit QED, is well described by the 
Hamiltonian 

H ujq X t , / X\ ] ^ t t (^\ 

— = —(Jz a^aaz-\-[uJc a'a aa (1) 

h 2 2 ^ \ 2) 2 ^ ^ 

taking into account only the lowest two energy levels of the 
qubit. The operators a^/a are the usual raising/ lowering 
operators for the harmonic oscillator and is the Pauli 
operator. In this description, we completely omit the mea- 
surement cavity as it is only used for reading out the state 
of the qubit and otherwise stays in its ground state. The 
energy level diagram described by the Hamiltonian given in 
Eqn. 1 can be seen in Fig. lb. The second term in Eqn. 1 
is the state-dependent shift per photon x/^tt = 9.4 MHz of 
the qubit transition frequency. The last two terms describe 
the cavity as an anharmonic oscillator with a dressed reso- 
nance frequency ujc and a nonlinearity K/27: = 325 kHz 
which is given by ~ x^/AKq [15]. All interaction 
strengths in the above Hamiltonian are at least one order of 
magnitude bigger than any decoherence rate in the system. 

To visualize and understand the evolution of the res- 
onator state, we measure the Husimi Q-function Qo in a 
space spanned by the expectation value of the dimension- 
less field quadratures Re{a} and Im{a}. Qo is defined as 
the modulus squared of the overlap of the resonator state 
1^) with a coherent state \a) hy QQ{a) = ^ \ {a\^) \^ . Al- 
ternatively, we can write Qo using the displacement oper- 
ator Do, = e"^^-"*^ (note that = D_c.) as Qo{a) = 
^ \ (0|I)_a|^) p which describes the actual measurement 
procedure employed in the experiment. The sequence to 
measure Qo can be seen in Fig. 2a. The initial displace- 
ment, 1^/3, creates a coherent state |^) = in the cavity, 
whose Qo is given by a Gaussian, ^e"'""^'^. After a vari- 
able waiting time t, we measure Qo{a) by displacing the 
cavity state by —a and determine the overlap of the re- 
sulting wavefunction with the cavity ground state. The 
population of the cavity ground state can be measured by 
applying a photon number state selective tt pulse, X^^^^ 
on the qubit (see supplementary material). The qubit is 
excited if and only if the cavity is in the n = Fock state 
after the analysis displacement. Applying tt pulses to the 
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FIG. 2. Technique for measuring the generalized Husimi 
Q-functions. a, The experimental pulse sequence consists of 
a 10 ns displacement pulse which creates a coherent state 
in the cavity. After a variable waiting time, t, we analyze the 
state in the cavity by displacing the cavity by —a followed by 
a TT pulse on the qubit conditioned on having n photons in the 
cavity. In this way we can measure the generalized Q-functions 
Qn(tt) = {n\D-a\^) 1^ which are projections of the displaced 
wavefunction onto the Fock states \n) where Qo{a) is the Husimi 
Q-function. b, Density plot of Qo of the ground state. We 
measure the Q-function at 441 different analysis displacements 
a. The tt pulse on the qubit X^^^ is conditioned on having no 
photons in the cavity after the analysis displacement, c, Linecut 
of Qo of the ground state along Im{a} = 0. The red line is a 
plot of the theory, a Gaussian given hy ^e~^^^ . A fit to the data 

with a Gauss function Ig-^^/^^^ results in 2cr^ = 1.03 ± 0.02 
which is consistent with the expected width, d. Density plot of 
Qi of the ground state. In this case the tt pulse on the qubit 
X^^-^ is conditioned on having one photon in the cavity after 
the analysis displacement, e, Linecut of Qi of the ground state 
along Im{a} = 0. The red line is a plot of the theory given by 
a Poisson distribution -laPe"'""' . 



qubit conditioned on other photon numbers, X^, allows 
us to measure the overlap of the displaced state with any 
Fock state n, which we will call the generalized Q-functions 
Qn{(^) = ^ \ In essence we can ask the ques- 

tion: "Are there n photons in the resonator?" , using photon 
number state selective pulses [26] . To test the analysis pro- 
tocol we measured Qo and Qi of the cavity in the ground 
state, Fig. 2b-e, by omitting the first displacement pulse of 
the sequence given in Fig. 2a. 

Using this method, we can follow the time evolution of 
a coherent state in the presence of the Kerr effect. In the 
experiment, we prepare a coherent state with an average 



photon number \/3\ 



n 



4 using a microwave pulse [23] 



to displace the cavity. We then measure Qo for different de- 
lays between the preparation and analysis pulses. A com- 
parison of the theoretical evolution of the coherent state 
and the measured evolution can be seen in Fig. 3. The 
time evolution of the state is described by considering the 
action of the Kerr Hamiltonian i^xerr on a coherent state 
1/3) in the cavity [5, 27]. In the rotating frame of the har- 
monic oscillator, with the qubit in the ground state, we can 



write 
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For short times, the nonlinear phase evolution of the Fock 
states |n) is closely approximated by a rotation of the state 
with an angle ^xerr = Kt{\l3\'^ + 1/2) with respect to the 
frame rotating at ujc. The onset of this rotation can be seen 
in Fig. 3a which is taken at the minimal waiting time of 15 
ns between the two displacement pulses. Due to this wait- 
ing time, the state rotates under the influence of the Kerr 
effect from /3 = 2 to /^e*"^^^"^^ = 2e*^--^^. For longer times we 
can see how the state rotates further and spreads out on a 
circle. Fig. 3b, c. This spreading can be simply understood 
in a semi-classical picture where the amplitude components 
in the coherent state further away from the origin evolve 
with a higher angular velocity given by the dependence 
of the Kerr effect. Complete phase collapse is reached at 
a time when the phase dispersion across the width of the 
photon number distribution corresponds to ~ tt, which can 
be estimated as Tcoi = [2]. For our system the com- 

plete phase collapse happens at 385 ns. Fig. 3b, shows 
that a Kerr medium can be used as a resource to generate 
squeezed states [27]. The state is squeezed along the Re{a} 
quadrature with a width of 0.88(1) as predicted from the- 
ory. The maximum squeezing occurs at t = 58 ns with a 
width of Qo of 0.87. 

After the complete phase collapse, structure re-emerges 
in the form of multi-component superpositions of coherent 
states. Fig. 3d-f, at times which are integer submultiples 
of the complete revival T^ev = Fig. 3h. The revivals, 
periodically appearing every T^ev, can be understood by 
noting that e*^"^ ^ = (— 1)"^ for t = T^ev The cavity state 
is then given by |^(T^ev)) = • At this time we get 
a complete state revival to a coherent state with opposite 
phase. For t = T^qv/q^ with q an integer larger than 1, we 
can write the state of the oscillator as a superposition of q 
coherent states [2] 



1 

2q 



2q-l 2q-l 
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For q=2 we get the two-component Schrodinger cat state 
similar to the cat states created in ref. [24, 28]. To dis- 
tinguish the q components of a cat state, the coherent 
states have to be separated by more than twice their width 
on a circle with a radius given by the initial displace- 
ment. In other words, the coherent states have to be quasi- 
orthogonal. This means for a displacement of = 2 the 
maximum number of coherent states that can be distin- 
guished is 4. 

In frame Fig. 3g we can see how the state again com- 
pletely dephases shortly before the coherent revival in 
Fig. 3h after t = 3065 ns. After this time we get a state 
with amplitude \f3\ = 1.78(2) which fits to the expected de- 
cay of the resonator state. The theoretical plots for Fig. 3 
were simulated by solving a master equation using the de- 
cay rate /^/27r = 10 kHz of the resonator and introducing 
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FIG. 3. Time evolution of Qo for a coherent state in the nonlinear cavity. Experiment (upper row) and theory (lower 
row) for a coherent state P — 2. The time t for the frames a-h is given above each panel. We measure Qo at 441 different analysis 
displacements a. The resolution of the pictures was doubled by interpolation. Initially the phase of the state spreads rapidly a-c, 
leading to a complete phase collapse after a characteristic time Tcoi = ~ 3^^^^- -^^^ short times, the Kerr interaction leads 

to a quadrature squeezed state along the Re{a} axis which can be seen in b. After the complete phase collapse structure emerges 
again, d-h, at times which are integer submultiples of the complete revival time Trev = At these times, one obtains coherent 
state superpositions which are multi-component cat states, up to a maximum number of resolvable components set by the average 
photon number in the initial coherent state displacement. The color scale of each pair of plots is individually rescaled with the 
scaling parameter given for a-h by A = 1, 0.93, 0.32, 0.24, 0.3, 0.46, 0.24, 0.89. The negative amplitudes in plot a are due to a 10% 
excited state population of the qubit which evolves in a different rotating frame (see supplementary material). This excited state 
population is not visible in the other frames as it disperses quickly and vanishes in the large positive amplitudes. 



a small detuning of 5 kHz of our drive from the resonator 
frequency ujc The time evolution of the state in the ex- 
periment agrees well with the theory. The hazy ring that 
can be seen in theory and experiment is produced by cav- 
ity decay. The evolution of the state from - 6.05 /is over 
50 frames, including two revivals, can be seen in a movie 
provided in the supplementary material. 

To get a more quantitative comparison of experiment and 
theory we want to determine the quantum state of the res- 
onator by measuring its Wigner function. Although in prin- 
ciple one can reconstruct the cavity wavefunction from the 
measured Husimi Q-function, in practice there is impor- 
tant information, such as the interference fringes between 
coherent state superpositions, which is exponentially sup- 
pressed as their separation increases [2] . This makes it hard 
to distinguish a mixture of coherent states from a coherent 
superposition. An experimentally more practical way to de- 
termine the quantum state of a resonator is to measure its 
Wigner function as it emphasizes the interference fringes. 
Previously, the Wigner function of a cavity [2] has been 
determined by measuring the parity of the resonator state. 
Here, we use a modified technique based on earlier work 
with ion traps and microwave circuits [3, 29]. The main 
difference is our ability to directly measure Qn{oi) which 
allows for a simple and efficient measurement and recon- 
struction of the density matrix of the resonator by using a 
least square fit to each Qn (see supplementary material). 
Using this density matrix we then calculate and plot the 
Wigner function. 

In Fig. 4 we show a comparison of the experimentally 
obtained Wigner functions to a simulation at three dif- 
ferent times during the state evolution. The times {t = 
27T/2K, 27T/3K, 271 /AK) were selected such that the Wigner 



functions correspond to two, three and four-component cat 
states. The simulation was again done by solving a mas- 
ter equation which includes the decay of the cavity. The 
fidelity F = (^id| Pm l^id) of the measured state pm, com- 
pared to an ideal n-component cat state |^id), consist- 
ing of coherent states with amplitude = 2e~^^l^^ is 
F2 = 0.71,^3 = 0.70,^4 = 0.71 for the two, three, four- 
component cats respectively. The Wigner functions show 
clear interference fringes which demonstrates that the evo- 
lution is indeed coherent and well described by the wave- 
function given in Eqn. 2, up to the decay of the cavity. The 
main reduction in the fidelity is due to the spurious excited 
state population of the qubit (see supplementary material) 
and the decay of the resonator. The decay of the resonator 
state is also responsible for the asymmetry in the interfer- 
ence fringes of the Wigner function, e.g. the maximum of 
the interference fringes for the two-component cat states is 
shifted to the left in both theory and experiment. 

In conclusion, we have shown that we can engineer strong 
photon-photon interactions in a cavity, entering the single- 
photon Kerr regime where K » We are able to ob- 
serve the collapse and revival of a coherent state due to 
the intensity-dependent dispersion between Fock states in 
the cavity. This opens the possibility to use such a Kerr 
medium for error correction schemes where a nonlinear cav- 
ity is used to realize the necessary components [9]. The 
good agreement between the theory and the experiment 
demonstrates the accurate understanding of this system. 
It also confirms our ability to predict higher-order cou- 
plings which is a necessary ingredient for understanding 
the behavior of large circuit QED systems. Furthermore, 
we have measured the evolution of a coherent state in a Kerr 
medium at the single photon level and shown a new exper- 
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FIG. 4. Wigner function of the multi-component cat 
states emerging during the Kerr interaction. The top 

row shows the measured Wigner functions of a coherent state 
subject to a Kerr interaction for a time t. The lower row shows 
the theoretically expected Wigner functions for the same in- 
teraction time obtained by a simulation including the decay 
of the cavity. The Wigner functions are reconstructed from 
measurements of the quasi probability distributions Qn(tt) for 
n = — 7. Each Qn{o:) was measured at the same displace- 
ments as the corresponding data in Fig. 3. Comparing the 
experimentally obtained state to an ideal cat state we find a 
fidelity F2 = 0.71,^3 = 0.70,^4 = 0.71 for the two, three, 
four-component cats respectively. The experimental data show 
excellent correspondence with theoretical predictions, including 
the interference fringes and regions of negative quasi-probability 
distribution, confirming that highly nonclassical states are pro- 
duced by the Kerr evolution. 

imental way for creating and measuring multi-component 
Schrodinger cat states. This demonstrates the ability to 
create, manipulate and visualize coherent states in a larger 
Hilbert space and opens up new directions for continuous 
variable quantum computation [30]. 
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CN 1 Materials and Methods 
> 

O 1.1 Qubit fabrication 
IZ 

Josephson junctions and dipole antennas are fabricated in one single fabrication step by electron-beam lithography 
followed by aluminum double-angle electron-beam evaporation. The two evaporations deposit thin aluminium films 
with a thickness of 20 nm and 60 nm respectively. These layers are separated by an AlOx barrier grown by thermal 
^ oxidation for 720 seconds in 2000 Pa static pressure of a gaseous mixture of 85% argon and 15% oxygen. 
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1.2 Measurement setup 

The state of the qubit is probed by measuring the qubit-state dependent cavity transmission using heterodyne detec- 



cd 

qH tion. Fig. 1 is a block diagram of the measurement setup. Microwave signals for control and measurement transmit 
through a coax cable with a 20 dB attenuator at 4 K and a 20 dB directional coupler followed by a 10 dB attenuator 
at 25 mK. The cavity transmission at the output goes through two cryogenic isolators at 25 mK followed by a HEMT 
^ amplifier with a noise temperature of 5 K at the 4 K stage. At room temperature, the transmission signal from the 
QQ cavity is further amplified by two low-noise room temperature amplifiers and is mixed down to a 20 MHz signal which 
fS| is digitized using a 1 GS analog to digital converter. 

(N 

^ 2 Full system Hamiltonian 

fSj The two cavities, one qubit system, can be described by the Hamiltonian 

> — = LJcCL^a + UJrnaL^m+(^qb^b (l) 

• ^ il 

in the strong dispersive limit of circuit QED. Where a/a^ , am/dlni ^/^^ raising/lowering operators for the storage 
cavity, the measurement cavity and the qubit. The Hamiltonian is calculated by treating the system as three coupled 
harmonic oscillators and introducing the cosine term in the Josephson relation, using the 0^ ~ (6 + 6^)^ term in the 
Taylor expansion, as a perturbation [1]. Higher order terms (~ 0^) are more than a factor of 1000 smaller for the 
parameters of this experiment and can be safely neglected. The Hamiltonian in Eqn.l of the main text can be recovered 
from the above Hamiltonian by neglecting all but the lowest two levels of the qubit (b^b^bb = 0) and replacing b^b with 
|(1 + cr^). Furthermore all terms ~ aj^a^ are zero as the readout resonator always stays in the ground state during 
the Kerr evolution. 

The first three terms in the above Hamiltonian describe each mode as a harmonic oscillator with lOc the resonance 
frequency for the storage cavity, cj^ the resonance frequency for the measurement cavity and cOq the ground to first 
excited state transition frequency for the qubit. The next three terms describe the anharmonicity of each mode, the 
anharmoncity of the storage resonator K, the anharmonicity of the measurement cavity and the anharmonicity of 
the qubit Kq. The last three terms describe the state dependent shifts of each mode due to the state of the other two. 
This means that the transition frequency of a resonator not only depends on the state of the qubit but also on the 
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Figure 1: Block diagram of the measurement setup. 



state of the other cavity. In our system ah of these state dependent shifts are much bigger then any decay rate. The 
state dependent shift of the storage cavity to the qubit is x, of the measurement cavity to the qubit is Xqm and of the 
storage cavity to the readout cavity Xcm- The values of ah parameters determined by spectroscopy and a comparison 
to the theory values obtained using finite-element calculations for the actual geometry, combined with "Black-Box" 
circuit quantization [1] are given in Table 1. 



Table 1: Comparison of the experimentally obtained and predicted values for the frequency and anharmonicity of each 
mode as well as the state dependent shifts between the modes. 





Exp. (MHz) 


Th. (MHz) 


deviation (%) 


(jUq /27r 


7850.3 


7890 


< 1 


cjc/27r 


9274.7 


9372 


1 




8256.4 


8336 


1 




73.4 


72 


2 


K/2ii 


0.325 


0.25 


30 




3.8 


3.7 


3 




9.4 


8.2 


15 


Xgm/27r 


29.5 


29.5 


< 1 


Xcm/27r 


2.45 


2.1 


16 



3 Qubit state readout and Cross Kerr dependence 

The state dependent shift between the two cavities, or Cross-Kerr, enables us to readout whether one cavity is in the 
ground state or not by probing the other cavity exactly like reading out the qubit state. This also means that the 
qubit cannot be readout independently of the state of the storage cavity. Both qubit and storage cavity in an excited 
state will contribute to the readout voltage. To remove this effect in our measurements, e.g. determining whether the 
qubit was excited by a photon number state selective tt pulse, we always perform a control experiment without a tt 
pulse applied to the qubit. This measurement gives us the readout voltage corresponding to the state of the storage 
cavity only. As our readout is linear, we can subtract this voltage from the combined measurement and thus infer the 
qubit state. This subtraction was done to measure all Qn{oL) in the main text. 
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4 Photon number state selective pulses 



In order to measure the state of the storage cavity, we must project the cavity onto a photon number (Fock) state. 
From Eqn. 1, we can see that the qubit transition frequency is dependent on the photon number state of the resonator. 
As long as the frequency shift, is much larger than spectral width a of the interrogation pulse, a photon number 
state selective pulse can be performed. In this experiment, x/27r = 9.4 MHz and cr/27r = 2.6 MHz, making a pulse 

selectivity of 1 — e 2<t2 > 0.99. This is essentially a CNOT [2] operation, with the qubit as the target, conditioned on 
the photon number in the cavity. This process can be described by 

oo ^ oo 

cavity) (8) {"^ qubit) ^^Pn W) ® \g) ^ Pn \n) \g) + Pm \m) ® \e) . (2) 

n=0 n=0 

Where the state of the cavity is described as a sum over the Fock states \n) with pn the complex amplitude. This 
procedure entangles the qubit and cavity thus a measurement of the qubit state will project the cavity onto the m*^ 
photon fock state. The probability of finding the qubit in the excited state is then given by Pm which is exactly the 
probability of finding the cavity in the state |m). 



5 Determining the cavity anharmonicity 

The cavity anharmonicity is measured by spectroscopy using a 50 ^s pulse to excite the cavity. The cavity excitation 
is measured by mapping the population left in the ground state to the qubit using a photon number state selective tt 
pulse. The tt pulse on the qubit will be successfull only when the cavity is in the ground state after the spectroscopy 
pulse. Spectroscopy on the cavity for a varying spectroscopy power is shown in Fig. 2a. For the lowest power the cavity 
is only excited on the |0) |1) transition at a frequency of 9.2747 GHz. As the power is increased different transitions 
appear in the spectroscopy. These transitions are multi-photon transitions from |0) to \n) with n = 2,3 which means 
the separation between the peaks is given by 0.5K/27r = 163 kHz, the Kerr nonlinearity of the cavity. The data were 
fit using a sum of Lorentzian peaks. 
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Figure 2: Cavity spectroscopy and displacement calibration, a. Cavity spectroscopy for varying spectroscopy 
power. The cavity excitation is measured by mapping the population left in the ground state to the qubit using a 
photon number state selective tt pulse. For the lowest power (o) the cavity is only excited on the |0) |1) transition. 
As the power is increased ( • ^ m ^ •) different transitions appear in the spectroscopy. These transitions 

are n-photon transition from |0) to \n) with n = 1,2,3 labeled from the right peak to left. The separation between 
the peaks is given by 0.5K/27r = 163 kHz. The dashed lines are fit to the data using a multipeak Lorentzian 
function, b. Photon number probability of the Fock states ri = . . . 7 (•, •, , ,•,•,•,•) for a 10 ns displacement 
with amplitude |a| of the cavity ground state. The axes were scaled by fitting a Poisson distribution with two free fit 
parameters to all seven Fock state populations simultaneously. The dashed hues are given by a Poisson distribution 
P(|a|) = |a|2^e-l"lV^! for the Fock state n. 



6 Cavity displacement & displacement calibration 

The cavity is displaced using a 10 ns pulse resonant with the |0) |1) transition of the cavity creating a photon number 
probability distribution given by a coherent state. We can measure the probability of finding the Fock state \n) by 
applying a photon number state selective tt pulse to the qubit and measure the qubit excitation. The displacement 
voltage e applied to the cavity is calibrated in units of the average photon number in the cavity n, or rather the 
corresponding displacement amplitude |a| = y^? by fitting a Poisson distribution P^(e) = (^)2^e('/^') /n!, with 
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the normalization constant Ae, to all measured photon number probabilities simultaneously. The probability of 
finding the Fock states n = ... 7 for a displacement with amplitude \a\ of the ground state can be seen in Fig.2b. 
The probabilities closely follow the expected Poisson distribution. This demonstrates that it is possible to create a 
coherent state despite the cavity non-linearity as long as the displacement pulse has a large enough spectral width. 

7 Excited state qubit population and negative values in the Hussimi 
Q-function. 

We found, that the qubit has a 10% excited state population which has been observed in other 3D superconducting 
qubits [3, 4] but is not yet completely understood. In our case this excited state qubit population leads to negative 
values in the Hussimi Q-function, as can be seen in Fig. 3a of the main text. These negative values are a result of the 
measurement procedure as Qo{c^) = is given by the readout voltage corresponding to the steady state population of 
the qubit. When the qubit excited state population is rotated to the ground state the readout voltage will be lower 
than for Qo{a) — 0. This rotation becomes possible as the resonator state for the qubit in the excited state evolves 
at a rate diff'erent by x compared to the ground state qubit resonator. The two coherent states will separate in phase 
space as can be seen in Fig. 3a of the main text. The evolution of this coherent state correlated with an excited state 
qubit will closely mirror the ground state population up to an additional rotation with x- This means that at T^e^/q 
we will also generated q-component cat states which have negative readout voltages. This will reduce the amplitudes 
and coherent fringes of the observed cat states as the amplitudes of our measurements are scaled to the amplitude 
of the coherent state of Fig. 3a in the main text. Effectively, the excited state population will lead to a reduction of 
the fidelity for the measured cat states. 

8 Cavity state tomography 

Using the photon number state selective pulse to project onto the n*^ photon Fock state, along with a preceding cavity 
displacement, we can measure (see Fig. 3): 

Qn{otm) = -{n\D{-am)pD{am)\n) (3) 

TT 

where n is the fock state projection; Q!^,the tomography displacement; and p, the cavity density matrix. The Qn are 
normalized probability distributions. Assuming a truncated photon basis, we are also able to directly reconstruct the 
cavity Wigner function [5, 6]: 

W{am) = - TT{D{-am)pD{am)P) 

TT 

= I (4) 

n 

= 2^(-l)"Q„(a„) 

n 

where P is the photon number parity operator. A more effective method for determining the Wigner function of the 
resonator is to reconstruct its density matrix. Realizing that the Qn are a set of linear equations, we can rewrite 
Eqn. 3 as: 

where Mnmij = {n\D{—am)\'i){j\D{c^m)\'n) and pij = {i\p\j). Measuring the distributions, Qn{c^m), with known 
photon projections n and displacements a^, we can perform a least squares regression to determine the cavity state, 
p. For this experiment, we measured 441 displacements each for projections onto the Fock states with to 7 photons, 
3528 measurements in total with 1000 averages per measurement point for each regression. We perform the linear 
regression with a priori assumptions that the cavity density matrix is Hermitian, of trace one, positive semi-definite 
and having a trunctated basis of ten photons. This reconstructed density matrix is then used to calculate state fidelities 
and plot Wigner functions. 
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Figure 3: Qnioc) probability distributions, a-h show the measured projections on the Fock states with n = — 7 for 
the two-component Schrodinger cat state. Each Qn{c^) was measured at 441 different displacements a. The different 
Qn{c^) functions were measured by conditioning the qubit tt pulse on the photon number in the cavity 
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